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We study the stability of static, spherically symmetric, traversable wormholes existing due to conformal continuations 
in a class of scalar-tensor theories with zero scalar field potential (so that Fisher's well-known scalar-vacuum solution 
holds in the Einstein conformal frame). Specific examples of such wormholes are those with nonminimally (e.g., 
conformally) coupled scalar fields. All boundary conditions for scalar and metric perturbations are taken into 
account. All such wormholes are shown to be unstable under spherically symmetric perturbations. The instability 
is proved analytically with the aid of the theory of self-adjoint operators in Hilbert space and is confirmed by a 
numerical computation. 



1. Introduction 

In our recent paper ^ we have considered spherically 
symmetric perturbations of wormhole solutions to the 
Einstein-massless scalar field equations which exist for 
scalar fields nonminimally coupled to gravity [21 El ■ The 
equations of motion were reduced to a single wave equa- 
tion for the scalar field perturbation which in this case 
comprises the only dynamical degree of freedom. An 
analysis of this wave equation leads to the conclusion 
that such wormholes are unstable under spherically sym- 
metric (monopolc) perturbations, and this instabilty is 
of catastrophic nature since the increment of perturba- 
tion growth has no upper bound. 

In this paper we continue this study and extend it in 
two respects. First, we discuss more general background 
configurations, namely, static, spherically symmetric 
wormholes which appear in arbitrary scalar-tensor theo- 
ries (STT) of gravity in which the effective gravitational 
constant can change its sign due to conformal continu- 
ation (CC) 4 . The investigation is, however, restricted 
to massless fields for which Fisher's well-known solu- 
tion holds in the Einstein frame. Second, we examine 
the problem in more detail, including the behaviour 
of metric perturbations related to those of the scalar 
field. A physically meaningful metric perturbation of 
an initially regular configuration should be regular ev- 
erywhere. This requirement turns out to impose an 
additional constraint on the scalar field perturbations, 
which makes the stability problem quite nontrivial. We 
finally prove that there exists at least a single growing 
mode of physically meaningful perturbations, i.e., such 
wormholes are indeed unstable. However, contrary to 
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the conclusion of Ref. pp, the perturbation grows at a 
finite rate. 

We thus find that gravitational instabilities, whose 
existence seems to be quite natural at surfaces where the 
gravitational coupling changes its sign (see, e.g., Ref. jSl 
for a discussion in a cosmological setting), still need 
much effort in their detailed study and even discovery. 

The paper is organized as follows. Sec. 2 is a brief 
description of the background static configuration and 
its place among more general configurations of this kind, 

1. e., static, spherically symmetric wormhole solutions of 
a general class of scalar-tensor theories (STT) admit- 
ting conformal continuations (CCs) Sec. 3 discusses 
spherically symmetric perturbation equations and the 
corresponding gauge freedom. Sec. 4 is devoted to a sta- 
bility investigation for wormholes, both analytical, using 
the theory of self-adjoint operators in Hilbert space, and 
numerical. 

2. Conformal continuations and 
wormhole solutions of scalar-tensor 
theories 

2.1. STT in Jordan and Einstein pictures 

Consider the general (Bergmann-Wagoner-Nordtvedt) 
class of STT, where gravity is characterized by the met- 
ric g^u and the scalar field cj); the action is 

s = j d^x^ifimig] + K4>)g'"'(|^,^.<ii^u 

-2[/(0) + 167rGi™}. (1) 

Here Tl[g] is the scalar curvature, /, h and U are cer- 
tain functions of (j) , varying from theory to theory, L„i 
is the matter Lagrangian, and G is the gravitational 
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constant, not necessarily coinciding with its Newtonian 
value. 

The action is simplified by the well-known con- 
formal mapping |H] 



(2) 



accompanied by the scalar field transformation (/) i—t ip 
such that 



In terms of y^^, and ip , the action for U = L.^ = , the 
case of massless scalar-vacuum fields to be considered 
here, takes the form 

[ d^x^{s\gnf){n[g] + [sign/(0)]r''V',^^.} (4) 



(up to a boundary term which does not affect the field 
equations). Here R\g] is the Ricci scalar obtained from 

The space-time Mj = M\g\ with the metric g^^, is 
referred to as the Jordan conformal frame (or picture) , 
generally regarded as the physical frame in STT; the 
Einstein conformal frame Me = with the field V' 
then plays an auxiliary role (see, however, discussions 
of the physical meaning of various conformal frames in 
13 IS] and references therein). The action corre- 
sponds to conventional general relativity if / > , and 
the normal sign of scalar kinetic energy is obtained for 
/((/)) > 0. Scalar fields in anomalous STT, in which 
< 0, lead to a kinetic term in Q with a "wrong" 
sign, are called phantom scalar fields. Such fields (with 
different potentials) are sometimes invoked in modern 
cosmological studies to describe dark energy. 

Exact static, spherically symmetric scalar-vacuum 
solutions of the theory ^ are well known Among 
them, wormhole solutions are generic in the case of 
phantom scalar fields [2]. Their stress-energy ten- 
sor T^i, manifestly violates the null energy condition 
(NEC) Tf,^k^'k'' > 0, kf^k^ = 0, such violation being a 
necessary condition for wormhole existence .10^ , there- 
fore wormhole solutions in their presence would have 
been naturally expected. One can note that, accord- 
ing to such solutions, both space-times Mj and Me 
have wormhole properties, i.e., represent regular static 
traversable bridges between two flat asymptotics. The 
stability of such configurations has also been proved 
[ 111 I12| by a direct study of perturbation equations, 
though it seems quite strange for a field system with 
energy density unbounded from below. This question 
evidently deserves further investigation. 

Our interest here will be in other wormhole solutions 
which appear in STT with > 0, in cases when the 
space-time manifold Me is mapped, according to (0), to 
only a part of Mj ; this phenomenon was named confor- 
mal continuation (CC) ^,13. The Jordan space-time 
Mj is then globally regular, represents a wormhole, and 



it two non-intersecting parts map to two singular space- 



times 



HE 



and 



HE 



2.2. Fisher's solution and its conformal 
continuations 

The general static, spherically symmetric solution to the 
Einstein-scalar equations that follow from Q), was first 
found by Fisher and was repeatedly rediscovered af- 
terwards. Let us write it in the form suggested in 
restricting ourselves to the "normal" case I > 0: 



^{u) = Cu + V'o, 



dt^ 



sinh (ku) 



k^du^ 
sinh^(fcu) 



dn-" 



(5) 



(6) 



where the subscript "E" stands for the Einstein frame; 
dfl^ — dO^ + siv? 9 d(p is the linear element on a unit 
sphere; C (scalar charge), m > (mass in geometric 
units). A: > and ipo are integration constants, of which 
the first three are related by 



k'=' 



(7) 



Without loss of generality we put C > and ipo=0- We 
are here using the harmonic radial coordinate u G in 
Me [g] , satisfying the coordinate condition a = 2f3 + j . 

Another convenient form of the solution is obtained 
in isotropic coordinates: with y = tanh(fcu/2) , Eqs. (0), 
(jnj are converted to 



'0(y) = 



c 



In 



1+2/ 



1-y 



dsl 



= A{y) dt^ 



V^A{y) 



[dy^ + y^dQ^ 



1 + y 



2m/ k 



(8) 



(9) 



The solution is asymptotically flat at u — > (y — > 
0), has no horizon when C ^ and is singular at the 
centre (u ~f oo , y— >1 — 0, ip —^ oo). When the 
scalar field is "switched off" (C = 0, k = m), the 
Schwarzschild solution is recovered. 

A feature of importance is the invariance of JSJ, lO 
under the inversion y i-^ noticed probably for the 
first time by Mitskievich T^. Due to this invariance, 
the solution 0, © considered in the range y > 1 de- 
scribes quite a similar configuration, but now y — > cx) 
is a flat asymptotic and y — s- 1 -f is a singular centre. 
An attempt to unify the two ranges of y , or, in other 
words, the two copies of Fisher's solution, is meaning- 
less due to the singularity at y ^ 1. We shall see that 
such a unification, leading to a wormhole, is achieved 
in Mj [g] where the singularity is smoothed out (in case 
C — \/6m) owing to the conformal factor. 

The corresponding Jordan-frame solution for any 
f{(p) and h{(j)) such that /(</>) > are obtained from 
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©, ijnji using 0, lO- If the function f{(j)) is everywhere 
finite, Mj has the same basic properties as Me . 

However, according to "1], there is a class of STT in 
which some solutions produce structures of Mj dras- 
tically different from that of Mg . Namely, let us use 
the field reparametrization freedom of the action 
[4> = 4>{4>ncw)] aud fix the paramctrization by putting in 
(CJ h{(j)) = 1. Then g], if the function /(0) has a sim- 
ple zero at some (j) — (pQ, there is a subfamily of static, 
spherically symmetric solutions to the field equations 
admitting a CC. The latter means that a singular sur- 
face in Me , corresponding to (j) ~ (j)Q , maps according 
to to a regular surface Strans in Mj . Then M can 
be continued in a regular manner through this surface, 
and the global properties of Mj can be considerably 
richer than those of Me : in the new region one can 
possibly find, e.g., a horizon or another spatial infinity. 
The above result was obtained in 4 for STT Q in ar- 
bitrary dimensions and with arbitrary potentials U{4)) . 
It was also shown 4 that a wormhole was a generic 
type of a conformally continued Jordan-frame manifold. 
Before studying perturbations of such generic solutions 
(but with U{4>) = 0, so that we have Fisher's solution 
in the Einstein picture), we discuss a specific example 
which makes evident the relations between Einstein and 
Jordan quantities. 

2.3. Example: wormholes with a conformally 
coupled scalar field 

A particular example of a CC is given by a free massless 
conformally coupled scalar field in GR. The latter is 
obtained when we put in 

/(<^) = 1 - 02/6, h{(j)) = 1. U{^) = L,„ = 0. (10) 

A transition sphere Strans, if any, corresponds to (j)^=Q. 
The transformation JJ)) now takes the form 



1-02/6 



(11) 



We assume that spatial infinity, where "0 ~* 0, cor- 
responds in the Jordan space-time Mj to < -\/6, 
where /(</>) > 0, so that the gravitational coupling has 
its normal sign. Then (|ll|l gives 

ip ^ Vg tanh"^((/)/\/6) + tpo, V'o = const. (12) 

Using ||SJ| and ©, it is now easy to write the metric 
in the Jordan picture. 

A CC through the sphere Strans (u = oo, y = 1, 
(f> ~ V6), which is singular in Me , is obtained when 
the infinity of the conformal factor 1// in compen- 
sates the zero of both and 'ggg simultaneously. This 
happens when, in accord with iQ, 



fc = 2m = 2C/V6, 



(13) 



which selects a special subfamily among all solutions. 
We will restrict the consideration to this subfamily. 



In terms of the isotropic coordinate y , the solution 
in the Jordan picture has the form [2] 



(1+TOo)' 



dt^ 



{dy^ + y^dQ^) 



1 + yyo 



(14) 



where yo = ta,nh{ipo / VQ) ■ The range < y < 1, de- 
scribing the whole manifold Me in the Fisher solution, 
corresponds to only a region Mj ' of the manifold Mj 
of the solution l|14(l . In all cases, y = corresponds to 
a fiat asymptotic, where </> V^yo < V6- The global 
properties of the solution depend on the sign of yo : 

yo < 0. The solution is defined in the range < 
y < l/|yo|- At y = l/|yo|, there is a naked attracting 
central singularity: ytt — > , — > , (p oo . 

b) yo = 0, (j) = V6y, y G M+ . In this case it is helpful 
to pass to the conventional radial coordinate r , substi- 
tuting y = m/{r — m) . The solution 



= (1 - m/rfdt^ 
(j) = V6'm/{r^m) 



dr^ 



(l-m/r)2 



(15) 



is the well-known BH with a conformal scalar field |15[ 
HQ). The infinite value of (f> at the horizon r = m does 
not make the metric singular since, as is easily verified, 
the energy-momentum tensor remains finite there. This 
solution has been shown to be unstable under radial 
perturbations |17|. 

c) yo > . This is the wormhole solution discussed in 
and, among other solutions, re-analyzed now. The 
solution is defined in the range y £ M.^ . At y ^ oo , 
we find another fiat spatial infinity, where (f> — > ^/6/yo , 
— ^ oo and gtt tends to a finite limit. 
The whole manifold Mj can be represented as the 
union 



IjlU Strans U] 



J1J2 



(16) 



where the region Mj ^ ( y < 1 ) is, according to Q , in 
one-to-one correspondence with the whole manifold Me 
of Fisher's solution jSJl, 10. The "antigravitational" 
(/(0) < 0) region Mj2 (y > 1) is in a similar cor- 
respondence with another "copy" of Fisher's solution, 
ME'[y], where, instead of (fT^ . 

= \/6 coth"^(</>/\/6) +i/'o, Vo^ const. (17) 

The metric g^^ of this second Einstein-frame manifold 
Me' should also be regularized by the factor 1/f on 
Strans , hcncc thc integration constants in it should sat- 
isfy the condition H13() . Moreover, one can verify that, to 
provide a smooth transition in the Jordan-frame metric 
y^i^ through Strans , all thc constants k , C and "00 
should coincide in Me and Me ' . The latter statement 
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is proved using the coordinate y which is common on 

both sides of Straus ■ 

This example weU iUustrates the general properties 
of conformally continued solutions ,4 . Namely, in the 
region beyond Strans i there can be a singularity due 
to = 0, as happens in the above case a) at y = 
I/I2/0I > 1- If there is no such singularity, we obtain a 
wormhole. Case b), with a horizon, is exceptional, in- 
herent only to the field (|10|l in four dimensions. Thus, 
for a more general action discussed in 1 , 18 , with 

f{<j))^l~^4>\ hi^) = l. = L,„ - 0, (18) 

with the coupling constant ^ > 0, in case ^ > 1/6 all 
solutions exhibiting a CC describe wormholes, whereas 
for ^ < 1/6 everything depends on an integration con- 
stant similar to yo in the above example, and we may 
have either a wormhole or a naked singularity. 

The stability analysis developed below covers worm- 
hole solutions obtained in the theory Q under the con- 
ditions 

h=l, U^Lm^O, Z((/')>0, (19) 

with an arbitrary function /(0), having a simple zero 
at some (f) = 4>o . In other words, we consider massless 
scalar fields in a general non-phantom STT, for which 
wormhole solutions exist due to a CC. 

3. Spherically symmetric perturbations 
and gauge freedom 

Consider small spherically symmetric (monopole) per- 
turbations of any static, spherically symmetric solution 
of the theory H19|l . The only dynamical degree of 
freedom is evidently related to the scalar field due to 
the generalized Birkhoff theorem jTH): if we take a time- 
independent scalar field, the equations of motion auto- 
matically lead to a static solution. 

We will use, for simplicity, the Einstein conformal 
frame, since the perturbation equations in Mj , being 
equivalent to those in Me , look much more complicated, 
and it is even hard to decouple different components of 
the Einstein equations. However, the boundary con- 
ditions that select physically meaningful perturbation 
should be formulated for variables specified in Mj and 
only then converted to Einstein-frame quantities. 

In the Einstein picture, the equations of motion are 
the Einstein-scalar field equations due to 

VaV"V = 0, (20) 

K = -t^t"- (21) 

We now write the metric in Me in the form 

dsl = e^'^dt^ ~ e^'^du^ ~ e^f^dn^, (22) 

where the functions a, /3, 7 as well as the scalar field 
if) are split into a static background part and a small 
(linear) time-dependent perturbation: 

a — a{u) + Sa{u, t) 



where u is a radial coordinate, and similarly for /3, 7 
and ■0 . Now, in addition to the freedom of choosing 
the radial coordinate u , we have an additional freedom 
of specifying the frame of reference in perturbed space- 
time, called gauge freedom. The latter makes it possible 
to specify (by hand) some linear relation between the 
perturbations. Certain care is needed to ensure that the 
resulting perturbation will not be a "pure gauge", i.e., 
will not be removable by coordinate transformations. 

Consider an infinitesimal coordinate transformation 
of the static metric ()22|l preserving its spherical symme- 
try, i.e., 

<ew = <ld + C^ = (^, C" 0, 0), (23) 

where the time dependence of the perturbations ^, ^ is 
separated: fj = ri{u) e^* and ^ = £,{u) e^^* . To preserve 
the diagonal form of the metric, we should put = 
^' g27-2a (where a and 7 are unperturbed), so that 
all perturbations are expressed in terms of rjiu) . For 
the metric functions and the scalar field we obtain 
(omitting the factor e^* ) 

Sa = ^e2T-2"[^" + ^'(2y -a')], 
Sj = l[f]2^-H7Ve2^-2"], 

S4> = ^V'e'^-'", (24) 

where the prime denotes d/du. 

Let there be a static configuration with (3' ^ and 
(f)' ^ . Then, if we have nontrivial time-dependent 
perturbations under the gauge condition 5/3 = (or 
Scf) = 0), Eqs. H24|l immediately lead to 77' = 0, which 
means that our perturbation cannot be caused by a 
transformation like 1)23(1 . i.e., is physical. The same 
is true for any gauge of the form fiSP + f2S(f> = 
where /i and /2 are any fixed functions of u, pro- 
vided fi(3' + f2<f>' 7^ 0. The reason is that (3 and (f> 
are scalars with respect to coordinate transformations 
of the 2-surfaces (a;°, x^). Thus, choosing such gauges, 
we can be sure that the perturbations to be studied will 
be physical. For other gauges, involving Sa and/or (^7, 
an additional inspection will be required. 

A more general approach to the problem of gauge 
in perturbation theory for spherically symmetric space- 
times can be found in Ref. ||2Qj; though, in the present 
case, our explicit treatment seems more transparent. 

4. Stability analysis 
4.1. The problem 

We have considered our set of linear perturbation equa- 
tions using two different systems of analytical computa- 
tion. Maple and Mathematica, which made it possible 
to compare the results and to avoid errors. 
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We use the Einstein conformal frame, in which the 
equations are much simpler, and the gauge 



(25) 



which is manifestly physical (see Sec. 3) and, in addition, 
transforms to Scj) — in the Jordan frame. Moreover, 
according to Eq. (|2U|I . we have the following relation be- 
tween the metric perturbations: 



5a = 2513 + Sj. 



(26) 



Two independent components of the Einstein equations 
for perturbations in the gauge (|25|) may be written as 

e^-'i?? = 2[Sl3' - P'i5p + S^) - j'SP] = 0, 
e'^"Rl = 2(3"{25(3 + 6j) 

- 2 e"'3+2^(5/3 + e'^I^Sp - 5p" = (27) 

Here primes denote derivatives with respect to u, the 
harmonic radial coordinate in the Einstein frame, a , (3 
and 7 describe the background configuration and satisfy 
the static field equations. We separate the variables 
using the substitution 



5{r,t) = 5{u) e^^*. 



(28) 



where (5 is a perturbation of any variable in our problem. 
After substitution of (|28|l into H27|) , 57 is expressed from 
the first equation, and then the second equation takes 
the form ( e^* is omitted) 

5/3" - 513 /(u) + F{u) 50 + G{u) 6(3 = 0, (29) 

where s, F, G are functions of u obtained from the 
metric 



F{u) = -2/37/3', 

G{u) = -2/3" + 2/3"77/3' + 2 6""^+^^, 
s{u) 



(30) 



A few words about the boundary conditions. At 
spatial infinity the choice is evident: 6(3 ^ 0. At the 
transition sphere Strans 6p should be finite, as well as 
its first two derivatives in u . This is necessary for the 
metric perturbations in the Jordan picture to be finite 
and smooth at Strans , which is easily checked using the 
transformation Q, @ and expressions in terms of the 
invariant length in the Jordan frame. 

As usual, we perform a transition from (|29|l to a 
Schrodinger-like form of the perturbation equation: 



d^y/dx^ + [E- V{x)]y{x 
where 

— / s'^du, 
m J 

/ 1 r 

Fdu 



0, 



(31) 



X 



5/3=fexp(-l 
V{x) = 2(/3,,//3,)2 - (/3,,, + 2/3,,7,)//3, 



+ 3/3,, + 5/3^ - 4/3,7, - 2m^ e 



2 „27-2/3 



(32) 



(33) 



where the subscript x denotes djdx and E = —m?QF'. 
The notations are chosen in such a way that the poten- 
tial V{x) and the "energy" E are dimensionless. The 
asymptotic forms of V{x) are 

V{x) K, 2/x'^ {x ^ CO, spatial asymptotic), 
V{x) « -1/(4x2) [x ^ 0, the sphere Strans). (34) 

Thus we have a quadratic potential well at Strans ; which 
is placed at x = by choosing the proper value of the 
arbitrary constant in the definition of x in Eq. H32() . 

The boundary condition at spatial infinity (u —^ 0, 
X — !■ oo) is ?/ ^ while the asymptotic form of any 
solution of (|31|l with < at large |x| is 



2/ ~ Cie 



mO|a;| 



Ci, 



const. 



(35) 



An admissible solution is the one with Ci = 0, with 
only a decaying exponential. 

At the other asymptotic, x — > 0, the condition that 
follows from the above continuity requirements reads 
y I \fx < oo whereas the solution of H31|l behaves as 



y w \/x{C3 + Gi Inx). 



(36) 



It follows that we must select the solution with C4 = . 
In other words, our problem is to find out whether there 
is a solution to the boundary-value problem for Eq. 131|l 
such that y as X 00 , J// ^/x < 00 as x — > and 
E < Q. In the remainder of the section we solve this 
problem. 

4.2. Summary of the solution 

We begin with a proof of the fact that the Hamilto- 
nian operator H related to Eq. 1)311) is self-adjoint and 
is bounded from below. To this end, we use an auxiliary 
operator T which has the same singularity at Strans 
as H . The one-sided boundedness indicates that the 
real part of the increment cannot be infinite. A fur- 
ther comparison of T and H shows that the continuous 
parts of their spectra coincide and lie in the non-negative 
part of the real number axis. So, if there are any solu- 
tions of our boundary value problem with i? < 0, they 
belong to a discrete spectrum. 

To prove the existence of a solution with i? < we 
use the well-known fact from quantum mechanics (its 
more general form is called the minimax principle) that 
the lower bound /io of the spectrum of an operator A 
is the infimum of the functional 



(V',AV'), 



(37) 



where the parentheses denote the scalar product (de- 
fined a bit later), the infimum is taken on the set of 
functions "0 which lay in the definition domain of T, 
and the norm — 1. Thus the value ('0,A'0) for 
any specified function -0 is an upper estimate for /^o : 
and if it is negative, then /Lto < 0. Functions which may 
closely resemble the unknown function that realizes the 
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infinium can give values of the functional (|37|l closest 
to . We guess such a function, which shows that the 
ground state of H lies below zero. This function is a 
ground state of a certain operator which is similar to H 
but simpler. 

4.3. The solution 

Consider the auxiliary differential equation 



1 



y{x) = Ey[x) 



(38) 



and investigate the question of self-adjointness of the 
Schrodinger operator 



Ty{x) 



dx"^ 



4x2 



y{x) 



(39) 



on the subset D{T) of real Hilbert space L2([0, oo)) 
such that, for y{x) e D{T) , (a) our boundary conditions 
(BCs) hold (so that, e.g., \y\l ^ < oo as x ^ 0) and 
(b) Ty{x) € L2 ■ The space i2([0, 00)) is a Hilbert space 
with an inner (scalar) product defined as the Lebesgue 
integral 



((/3,'0) = / (p*^dx, 



(40) 



where the star stands for complex conjugation. D(T) is 
dense in L2 since C^(0,oo) C D(T), where C^(0,oo) 
is the the subset of functions in C°°(0, 00) with a com- 
pact support separated from 0. It is a dense subset in 
L2 EH. 

One can show that the operator defined in this 
way is symmetric and therefore closable Obviously, 
the BCs of our Hilbert space are homogenous. The 
Schrodinger equation (|38|l related to the operator T has 
the solution 



Ciy/xKo{V—Ex) + C2y/xIoiV—Ex), 



(41) 



where E is the "energy" corresponding to —m^il^ of 
our problem, so to prove the instability we should show 
that there are "quantum states" with E < 0; Kq and 
Iq are the zero-order modified Bcssel functions of the 
first kind. Neither of these functions, nor their any lin- 
ear combination, satisfy our BCs. This means that the 
operator T — EI , E < has a bounded inverse op- 
erator (T — EI)~^ with a definition domain dense in 
L2 ■ The existence of a reverse operator follows from 
the well-known alternative: under given homogenous 
BCs, either the differential equation L[y] — g{x) has 
a uniquely defined solution y{x) , or the homogeneous 
equation L[y] =0 has a non-zero solution. In our case. 



1 



(42) 



The boundedness and the density property of the defini- 
tion domain of [T — EI)~^ in L2 follow from studying 
the properties of solutions to the equation L[y\ — g{x) 



with nonzero g{x) £ L2 ■ The existence of (T — EI)^^ , 
E < Q means that the domain (—00, 0) C p{T) , p(T) 
being the resolvent set of T . 

Considering in a similar way Eq. (|42() with _E > , 
one can show that [0,oo) is a continuous spectrum. 

Thus we have shown that T is a closed symmetric 
operator which contains real numbers in its resolvent 
set. It satisfies the conditions of the second corollary 
of Theorem X.l in 121!: // the resolvent set of a closed 
symmetric operator contains at least one real number, 
then this operator is self-adjoint. The self-adjointness of 
this operator was also mentioned in passing in Ref . . 

The proved properties of T make it possible to use 
the wealth of results obtained in the theory of self- 
adjoint operators. In particular, we use the following 
two theorems: 

Theorem 1 (Rellich j^). Let A be a self-adjoint op- 
erator on D(A) and B a symmetric operator on D{B) , 
so that D{B) D D{A) and 

WBtPW < a\\iP\\+b\\Ai^\\, 

b < 1. Then the operator A + B is self-adjoint and 
D{A + B) = D{A). 

Theorem 2 (Kato [221 ). Let the conditions of The- 
orem 1 hold and A be bounded from below (or from 
above, or from both sides), then A-\-B is bounded from 
below (or from above, or from both sides), but not nec- 
essarily with the same bound (bounds). 

Considering T as A in these theorems, we can 
rewrite H31|l as 



Ty{x) + Vix)yix) = Eyix), 
where 

V{x) = V(x) + l/4x^ 



(43) 



(44) 



Since V is bounded (this is true since V{x) — s- as x — s- 
00, and V{x) is bounded everywhere), the conditions of 
Theorem 1 are fulfilled, and the operator 



H 



-d^y/dx^ + V{x)y{x) 



(45) 



connected with Eq. H31|l is self-adjoint on D{T) . 

Using the spectral theorem for unbounded operators 
|21| . we prove that T is non- negative and consequently 
is bounded from below. Therefore the operator H is 
bounded from below too (Theorem 2). 

We now wish to show that the continuous spectrum 
of H coincides with the continuous spectrum of T . We 
use the following theorem: 

Theorem 3 PI]. Let A he a self-adjoint operator 
and C a symmetric operator such that C{A'^ -I- i)^^, 
n S N is a compact operator. Then, if B = A + C is 
self-adjoint on D{A), cross(^) = cross(S)-^ 

^We denote: (Tdisc = discrete spectrum, (Tcont = continuous 
spectrum, (Tcss = essential spectrum 1211 . I n our CcisG, Ccss con- 
sists of cTcont and a possible limiting point of (Tdisc ■ 



Conformal continuations and wormhole instability in scalar-tensor gravity 



7 



The compactness of the operator V{T + i) ^ follows 
from its integral representation: 



{ViT + ^)-'f){x) 



dy K{x,y) f{y), 



K{x,y)^V{x)G{x,y), 



(46) 



where G{x, y) is the Green function, core of the integral 
operator {T + . As follows from the asymptotic 
properties of V{x) and G{x,y), K{x,y) £ L2([0,cxd) x 
[0,00)). So, V{T + i)~^ is a Hilbert-Schmidt operator 
and hence is compact [2] ■ The conditions of Theorem 
3 arc fulfilled, and (JassiH) — acssiT) = [0,cxd). Since 
— fcssUcTdisc EBi remaining part of cr{H) belongs 
to the discrete spectrum. 

Eq. (|29|l may also be expressed in a non-Schrodinger 
form using the new variable w instead of u : 



1 



(47) 



and further converted to a normal form: 

11 + 32^; 64 1 
1-T + 



le^; 



2w 



1 



3(1 + 4u;) 
32 



4(1 



6(1 + it;) (l + 4w)2 
where 

n{w) = S(3{w) exp 



n(w) + 



{2w - 1) dw 



-0, (48) 



(49) 



2it;(l + w){l + Aw)^ 
We cannot solve this equation, but its truncated version 
d'^n{w) . ll + 32£; 1 



has the solution: 

n{w) = ciM 



2w 



(11 + 32£;) 



4^2 



n{w) = (50) 



C2W 



(ll + 32£;) 



(51) 



where M and W are Whittaker functions. If the first 
argument, usually denoted as the index fc, takes on the 
values k = 1/2, 3/2, . . . , n+ 1/2, n £ {0,N}, then 
M{k, 0, w) = W{k, 0, w), and the corresponding quan- 
tity E has the values 



11 fc2 

^^ = "32~y'^ 



1 



11 

2fc2 



which satisfy the equation 
(11 + 32Ek) 



= k. 



(52) 



(53) 



Performing the inverse coordinate transformation 
from the variables of (|48|l to the variables of (|31|) , we see 
that the functions y corresponding to W{k, 0, Sy^—E/^w) 
belong to D{T) . 



Hence E takes a discrete set of values in the in- 
terval [—3/8 + \/23/32,0), and they lie in ddisc of the 
Schrodinger operator related to Eq. (|50|l . The latter op- 
erator has an infinite discrete spectrum with the limiting 
point . Any other eigenvalue E gives a solution related 
to the essential spectrum, or the resolvent set of the lat- 
ter operator. The solutions obtained make it possible 
to apply the minimax principle ( 21 , theorem XIII. 1). 
The theorem is applicable because H is self- adjoint and 
bounded from below. 

We use the solution with k = 1/2 as a trial function 
to find an upper estimate of the lower bound of a{H) : 



(54) 



where the parentheses denote the scalar product on L2 , 
i.e., Lebesgue integration, ip(Ei/2, x) is a normalized 
function obtained from W{l/2, 0, 8^ 
substitution which transforms Eq. 



-Ei/2w) by the 
into (EU. Ex- 
plicit integration in x is impossible because we cannot 
represent w{x) in elementary functions, but we can in- 
tegrate in w after necessary substitutions. It is con- 
venient to represent H as a sum of two Schrodinger 
operators, where the first, H' , corresponds to Eq. H5U|I 
and the second, V' , corresponds to the remaining part 
of the operator. We obtain 

^i'o^{ml/2,x),H' + V' ^{E,/2,X)) 
= E,/2 + mE,/2,x),V' i^{E^/2,x)) 

JW{ll2,Q,^^/^E^2wfU{w)dw 

= Ei/2 



J W{l/2,0,8^f^WJ^wyr{w)dw 



(55) 



where 

U{w) 



64 



3(1 + Aw) 



4(1 
1 



32 



^ 6(l + w) ^ (H-4u;)2' ^^^^ 

and r{w) — 16(l-|-l/w). Integration gives /ij, ~ —0.039. 
Since the essential spectrum begins from 0, according 
to the minimax principle, /Iq is an upper estimate of the 
ground state which is thus located below zero. 

Thus we have proved the existence of negative eigen- 
values of the operator H under physically justified BCs, 
and therefore there are exponentially growing solutions 
(at least one) of our perturbation equation (|29|l . 

We have also solved our boundary-value problem 
1)29(1 numerically, applying the Fortran procedure SLEIG 
and using the coordinate transformation (|47|) . We found 
a single discrete eigenvalue at —0.048, which fits our es- 
timate of /ip . The corresponding problem for Eq. I|48(l is 
not suitable for using SLEIG because Fortran docs not 
"understand" the BC y/\/x < 00. 

Recall that we have been working in the Einstein 
frame, so that the coordinates used cover only half of the 
wormhole, so that we should use two copies of this patch 
and verify whether the metric perturbations remain re- 
ally smooth at the transition sphere §tians ■ Compu- 
tation of the corresponding metric perturbations ((5/3, 
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57) at the point where /((/)) — shows that their first 
derivatives in I (where I is the Gaussian radial coordi- 
nate in the Jordan picture, such that gu = —1) are 
zero, while the second derivatives are finite, so that 
the linearized gravity equations are there meaningful 
and hold. We conclude that the metric perturbations 
found are physical, and the wormholcs under consid- 
eration are unstable under small spherically symmetric 
perturbations. The decay rate depends on the value 
of TO in ©, since E = —rri^^l^. The wormhole ra- 
dius (understood, for simplicity, as ^/— 322 at the tran- 
sition sphere rather than throat radius which is smaller 
but generically of the same order) is proportional to to 
and also depends on the integration constant in the so- 
lution of (PJ. Let us discuss the special case of con- 
formal coupling. According to p4|l . the wormhole ra- 
dius is = 2to-\/ (1 + ?;o)/(l — yo) ■ If we assume 
2/0 1 , then both rwh and the throat radius are ap- 
proximately equal to 2m. The characteristic time of 
decay, r = , is proportional to to (which has the 
dimension of length) : 



T ~ to/y^ ~ 5m. 



(57) 



For a wormhole radius of the order of a typical stellar 
size ~ 10^ km, the time t is a few seconds (slightly 
greater than the time needed for a light signal to cover 
the stellar diameter). If yo increases under fixed worm- 
hole radius, then m decreases, so r decreases too. We 
see that, for all wormholes with a fixed radius, r < 5m. 

Similar estimates can be obtained for other STT 
characterized by different f{(p) . 
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